A NOTE ON THE PAULI PROBLEM IN LIGHT OF 
APPROXIMATE JOINT MEASUREMENTS 



JUSSI SCHULTZ 

Abstract. We show that there exist informationally incomplete 
phase space observables such that the Cartesian margins are infor- 
mationally equivalent with position and momentum. This shows 
that it is possible to reconstruct the position and momentum dis- 
tributions of a quantum system from the statistics of a single ob- 
servable, and thus a single measurement, even though the state of 
the system is not uniquely determined by the statistics. 



I. Introduction 

The Pauli problem is the classic example concerning the possibility 
of performing quantum state tomography. In a footnote in [I], Pauli 
noted that the problem whether or not the state of a quantum system is 
uniquely determined by the position and momentum distributions "has 
still not been investigated in all its generality". In modern terminology, 
this is a question about the informational completeness |2] of the pair 
(Q, P) of position and momentum observables. It is now known that 
the answer to this question is in the negative: (Q, P) is informationally 
incomplete. As a matter of fact, a wide variety of counterexamples 
have been constructed showing that different states can have the same 
position and momentum distributions [2j [3j HI [5] . This has been viewed 
as "an illustration of the surplus information J6] coded in a quantum 
(pure) state when compared with its classical counterpart" [7J. 

The purpose of this paper is to address the Pauli problem in a slightly 
stronger form. We pose the following question: 

(Q) If the position and momentum distributions of a quantum sys- 
tem can be reconstructed from the statistics of a single observ- 
able, is the observable necessarily informationally complete? 
We show that the answer to (Q) is also in the negative. We con- 
struct explicitly a phase space observable whose Cartesian margins are 
informationally equivalent [8j with position and momentum, thus al- 
lowing the reconstruction of the corresponding distributions from the 
marginal statistics, even though the observable is informationally in- 
complete. This result further illustrates the difference in the role of 
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position and momentum in determining the state of a system in quan- 
tum mechanics, as opposed to that in classical mechanics. Indeed, the 
state of a classical system is given by its generalized position and mo- 
mentum coordinates which can, in principle, be determined by a single 
measurement, whereas in quantum mechanics even the simultaneous 
determination (in the sense of (Q)) of position and momentum does 
not guarantee unique state determination. 

The paper is organized as follows. In Section [2] we lay out the general 
framework for this study. We review the relevant definitions and results 
concerning informational equivalence, informational completeness and 
phase space observables. In Section [3] we give the counterexample to 
(Q). We also prove the existence of informationally complete phase 
space observables whose margins do not suffice to determine position 
and momentum. In Section H] we consider the process of reconstructing 
position and momentum distributions from the statistics of a single 
phase space observable. We do this using two different methods: the 
Fourier theory and the method of moments. The conclusions are given 
in Section 



2. Preliminaries 

Let H = L 2 (M.) be the Hilbert space associated with a quantum 
system such as a spinless particle confined to move in one dimension 
or a single mode electromagnetic field. Let {h n \n — 0, 1, 2, . . .} be the 
orthonormal basis of H consisting of the Hermite functions and denote 
V = \m{h n \n = 0, 1, 2, . . .} so that V = U. 

The states of the system are represented by positive trace class oper- 
ators p with unit trace. The observables are represented by normalized 
positive operator measures E : B(M. n ) — > C(H) where B(M. n ) denotes 
the a- algebra of Borel subsets of M. n , and C{%) is the set of bounded 
operators on H. An observable is called sharp if it is projection valued, 
that is, E(X) 2 = E(X) for all X e B(R n ). For a system in a state p 
the measurement outcome statistics of an observable E is given by the 
probability measure : B(W n ) -»■ [0, 1], p E p (X) = tr[pE(X)). 

Definition 1. Let E : B(R n ) C{H) and F : B(R m ) -»■ C(H) be 
observables. 

(a) If for any two states p and a, p^ = p^ implies p F p = p^, then 
the state distinction power of E is greater than or equal to that 
of F. 

(b) If for any two states p and a, p^ = p^ if and only if p^ = p^ ; 
then E and F are informationally equivalent. 
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(c) If for any two states p and a , p E = p E implies p = a, then E is 
informationally complete. 



Let Q and P denote the selfadjoint position and momentum opera- 
tors and let Q, P : B(R) — > Ci^H) be the corresponding sharp observ- 
ables. Define the Weyl operators W(q,p) = e l ~ e~ iqP e ip ® , (q,p) G M 2 , 
and for each positive trace class operator T with unit trace define the 
phase space observable G T : £>(M 2 ) — > C(H) by 

G T (Z) = ±- [ W(q,p)TW(q, P ydqdp 

for all Z G B(M. 2 ). The operator T is called the generating operator 
of G T . The Cartesian margins of G T are the unsharp position and 
momentum observables pF * Q, v T * P : B(M.) — > C(H) defined as 

(p T *Q)(X)= [ p T (X -x)dQ(x), 
Jx 

( 1 y T *P)(Y)= [ u T (Y-y)dP(y), 
Jy 

for all X, Y G B(W). The convolving measures are determined by the 
generating operator T so that p T {X) = tr[TQ(— X)] and v T (Y) = 
tr[TP(-y)]. 

In general the state distinction power of, say, the unsharp position 
p T *Q does not exceed that of Q. However, they may be informationally 
equivalent. This is the case if and only if the support of the Fourier 
transform of the convolving measure is R [9j, that is, supp/i T = M 
where 

(1) p T {p) = -L / e~^dp T {x) = 4=tr[TW(0,p)] 



for all j) 6 1. The same is of course true for the unsharp momentum 
v T * P, where the Fourier transform is now given by 

(2) u T (q) = -^tx[TW(-qM- 



The informational completeness of the phase space observable G T 
can also be characterized in terms of the Weyl transform (q,p) (->■ 
tr [TW(q,p)} of the generating operator. Indeed, a sufficient condition 
for informational completeness has been known since [10]. In the recent 
paper this question was exhaustively resolved and several equiva- 
lent necessary and sufficient conditions were obtained. For the purpose 
of the present paper the following characterization is convenient (see 
[IH Prop. 4]). 



Theorem 1. The phase space observable G T is informationally com- 
plete if and only if the support of (q,p) h-> tr [TW(q,p)] is M 2 . 

An important class of informationally complete phase space observables 
are those whose generating operator satisfies tr [TW(?,p)] 7^ for all 
(q,p) € IR 2 - Such operators are called regular |12] . Clearly any T 
determined by a Gaussian wavefunction is regular. 



3. Main results 

The phase space observables G T are archetypes of approximate joint 
observables for position and momentum. It is even possible to choose 
G T in such a way that that the Cartesian margins fi T * Q and v T * P are 
informationally equivalent with Q and P. Indeed, consider the simplest 
case T = \h )(h \ so that for any state the corresponding probability 
density is the Husimi Q-function [T3] of the state. In this case the con- 
volving measures are given by fi T (X) = v T (X) — f x e~ x dx so that 

W(v) = v T (p) = -^^e * which confirms informational equivalence. 
However, this particular observable is also informationally complete, 
as can be seen from (ho\W(q,p)h ) = e~^ q +p The following Propo- 
sition shows that this is not in general the case. 

Proposition 1. There exist informationally incomplete phase space 
observables whose margins are informationally equivalent with position 
and momentum. 

Proof. Let (p = Xi-1/2,1/2}, the characteristic function of the interval 
[— |], and define 

T=\\<p)(<p\+\mm 

where (p denotes the Fourier transform of (p, that is, 

ftp) = -L fe*<p{x)6x 1 Sin(p/2) 



for all p 6 M \ {0} and <p(0) = -h=. Now consider the observable G 5 



^tt J V2n P/ 2 

-T 

siuei nit; uusei vciuit: v. 

and, in particular, the margins fi T * Q and v T * P. Using the fact that 
ip and tp are even functions we obtain the convolving measures 

/i T (X) = u T (X) = tr[TQ(X)] 

i(|v.(x)| 2 + |^)| 2 ) dx 
x / 
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for all X G B(W). Equations ([T]) and ([2]) now give the Fourier trans- 
forms as (see Figure [T]) 



(p) = v (p) 



1 



2V27T 



e ipx \<f{x)\ 2 dx 



2V27T 



<^(:r)93(x + p) dx 



2 ^(l-H + ^) ; whenH<l 

1 sin(p/2) 
2V2^ P/2 ' 



otherwise. 



This shows that Ji T (p) = if and only if p = 2mr, nGZ \ {0}, that is, 



supp /i T = supp z/ T 




FIGURE 1 . The plot of J1 T as a function of p. 



Now consider the informational completeness of G T . Using the fact 
that F~ 1 W(q,p)F = W(-p, q) for all (q,p) G M 2 , where F is the 
Fourier- Plancherel operator, we get 

tr[TW(q,p)\ = \e- 1 ^ [ e ipx ^(x)(p(x - q) dx 



, 1 i^E 

+ 2 62 



y(x)<^(x + p) dx. 



The first integral is zero when \q\ > 1 and the second one vanishes for 
\p\ > 1. Thus ti[TW(q,p)} = whenever \q\, \p\ > 1 which shows that 
the support of (q,p) h-> tr[7W(g,p)] is not the whole M 2 . According to 
Theorem [T] this means that G T is not informationally complete. □ 



Note that since the generating operator T is typically linked to the 
state of the probe system used in the measurement of G T , the generating 
operator of Proposition [1] is also physically meaningful. Indeed, in the 
case that T is (mathematically) the state of the probe, this can be 
realized by randomly preparing the probe system in each experimental 
run into the state \(f){<p\ or \(p)((p\. That is, in each run the probe state 
is localized either in position or momentum space with probability 1/2. 

We close this section by showing that the informational completeness 
of a phase space observable does not imply that the margins are infor- 
mationally equivalent with position and momentum. In this case it 
seems difficult to give a physically relevant counterexample, but math- 
ematically it can be constructed. 

Proposition 2. There exist informationally complete phase space ob- 
servables whose margins are not informationally equivalent with posi- 
tion and momentum. 

Proof. To begin with, we pick any regular generating operator To, i.e., 
one that satisfies tr [T W(q,p)} ^ for all (q,p) G M 2 . Then for any 
positive / G L 1 (R 2 ) such that J f(q,p) dqdp = 1, the trace class oper- 
ator 

f*T = J f(q, P )W(q,p)T W(q, P ydqdp 

is positive and of unit trace [12] . Furthermore, we have tr [/ * T W(q, p)} 
2n f(—p,q)tr [ToW(q,p)] so that in view of Equations (TjQ) and (T5]), 
and Theorem [1] we need to find the function / in such a way that 
supp/ = M 2 but the supports of the functions p t— > f(—p,0) and 
q i — y f (0, — q) are proper subsets of IR. 

To begin the construction, we define for any r > the function 

9o,r^P) = ^- r e- q2 (x[- ilS ]*X[- Sli ])(p) 

so that 

1 1 — cos(ry) 

9oA^,y) = ^^ — - 2 — e 4 - 

In particular, g 0tr G L^IR 2 ), g , r >0,f g 0>r (x,y) dxdy = 1 and^ ,r(g,p) = 
if and only if \p\ > r. Then for any 9 G [0, |] define 

9e, r {x,y) = go tr (xcos8 + ysm9, —xsinO + ycos9), 

so that 7)q t is nonzero only on the strip which is in an angle of 9 with 
respect to the g-axis (see Figure [2]). In particular, ^e»,r 0) = when 
\q\ > ±g and g e , r { ,P) = when \p\ > ^e- 



FIGURE 2. The shaded region represents the strip in 
phase space where the function gQ r is nonzero. 



We then define 6„ = 



n 



and r„ = 



sin 9 n 
2 n -l 



and the function g via the 



series 



oo 2 n -l 



1 



g(x,y) = ^2Yl 



9ke n ,r n (q,p) 



n=l k=l 



so that g G L^IR 2 ) and g(x,y) = if and only if gke„,r„(x,y) = 
for all n G N and k = l,...,2 n — 1. Finally, we define f(x,y) = 



L^M 2 ), f(x,y) > for all (x,y) G 1R 2 , and / f(x,y) dxdy = I. 

In order to see that this function has the desired properties, let B be 
a ball with radius e > centered at (<7o,Po) £ Due to the symmetry 
properties of / we may without loss of generality assume that (fchPo > 0. 
Since {k6 n | n G N, k = 1, . . . , 2 n — 1} is dense in [0, |], there exists an 
n G N and a k G {1, . . . , 2™° — 1} such that the line p = gtan(/c 6' no ) 
intersects with B. In particular, (r cos(/co6'no)? r o sin(/c ^n )) £ -B for 
some r > and since /(r cos(A; 6'„ ), r sin(A;o6' no )) ^ 0, we conclude 
that supp / = 1R 2 . However, f(q, 0) = when 



C(g(x,y) +g(-x,y)), where C 



\ (fg{x,y) dxdy) 



so that / G 
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and f(0,p) = when 

\p\ > sup ^ = sup (^^^ = -, 

n,k n,k v / & 

so that the support of neither of these "marginal" functions is R. Hence, 
the observable G^* T ° is informationally complete but the marginal ob- 
servables are not informationally equivalent with position and momen- 
tum. □ 

Remark 1. Obviously the counterexample for Proposition [1] could also 
have been given via a similar construction. Indeed, if To is as before and 
we define f = C$ (go,i + 9n/2,i) where Co > is the normalization con- 
stant, we find that G^°* T ° is informationally incomplete but the margins 
satisfy the appropriate conditions for informational equivalence. How- 
ever, we found it useful to present the given counterexample since it is 
physically more relevant. 

4. Reconstructing position and momentum 

We will next demonstrate two methods for determining the position 
and momentum distributions from the marginal statistics of a phase 
space observable. The first one uses the theory of Fourier transforms 
and second one uses the statistical method of moments. For the first 
part only the minimal requirement of informational equivalence for the 
margins is required as for the second part we need a stronger condition, 
namely, that of exponential boundedness. 

4.1. Fourier theory. Let G T be such that p T * Q and v T * P are 
informationally equivalent with the corresponding sharp observables Q 
and P, that is, supp p T = supp v T = R. Let p be an arbitrary state and 
consider the first marginal observable. By taking a Fourier transform 

of the probability measure we obtain (p T * pS)(p) = Vz7rju T (p)p^ (p) 
for all pel. It follows that for all p G R such that p T {p) ^ we have 

-Q/ N 1 (^*P?)(P) 

(3) P?(p) : 



Since we know that the Fourier transform is a bounded continuous 
function, we can take appropriate limits on the right-hand side of Eq. 
([3]) to determine for those p G R for which p T (p) = 0. Since the 
Fourier transform is injective, this uniquely determines the measure p^. 

If we want to invert the Fourier transform to get explicitly the form 
of the position distribution, we need to assume that the right-hand side 
of Eq. ([3]) is integrable. To that end, suppose that p is a pure state 



given by a unit vector ip £ T>. The position distribution x t— > \ip(x)\ 2 
then belongs to the Schwartz space <S(R), and hence also the Fourier 
transform is in <S(R). Furthermore, since S(R) C L X (R) we may 
use the Fourier inversion theorem to obtain the position distribution 

,2 _J_ [ „,J.f*P*)(v) 



for almost all x e R. A similar treatment shows that we can reconstruct 
the momentum distribution from u T * p^. The position and momentum 
distributions can therefore be reconstructed explicitly at least for the 
dense set of vector states T>. Furthermore, this method works also 
for any finite mixture p = 5^_ o c n |'0n)('0n| where (ip n )t=o c ^ an d 
(c n )^ =0 C [0, 1], Xln=o c « = 1- ■'■ n ^ ner words, for any state p whose 
matrix representation with respect to the basis {h n \n = 0,1, 2,...} is 
finite. It is known that such states are dense in the set of all states 
(this follows, e.g., from 1 1 1. Theorem 1]). 

4.2. Method of moments. The statistical method of moments was 
presented in this context in [15] and was further illustrated in |16j . 
The idea is to reconstruct the moments of the position and momen- 
tum distributions from the moments of the marginal statistics. In this 
case the informational equivalence of the observables is not sufficient 
to ensure the existence of the moments. Indeed, in the counterexample 
of Proposition [U the convolving measures p T and v T do not have any 
finite moments. Moreover, even the existence of finite moments does 
not guarantee that the probability measure is uniquely determined by 
the moments. In this sense the ability to reconstruct the moments of 
the position and momentum distributions does not necessarily mean 
that the actual distributions can be reconstructed. To circumvent this 
problem, we need to make a stronger assumption of exponential bound- 
edness for the measures. 

Recall that a probability measure /i : S(R) — > [0, 1] is exponentially 
bounded if there exists an a > such that 



According to [TTf Prop. 2] (which is based on the proof of [TS| Prop. 
1.5]) a probability measure is exponentially bounded if and only if there 
exist positive constants C, R > such that 

(4) | n[k] | < CR k k\ 

for all fceN. In addition to the existence of all moments, exponentially 
bounded measures have the important property of being determinate. 



That is, if \i is exponentially bounded and v : £>(R) — > [0, 1] is another 
probability measure such that u[k] = p[k] for all k G N, then v = p. It 
is also worth noting that any exponentially bounded measure satisfies 
supp/i = R. This is due to the fact that the Fourier transform has 
an analytic continuation to the strip {z G C | |Im(z)| < a}, and is 
therefore zero in at most countably many points. In particular, if p T 
and v T are exponentially bounded, then the corresponding observables 
p T * Q and v T * P are informationally equivalent with Q and P. 

Let again p be a state and suppose that G T is such that the mar- 
ginal measures p T * and v T * p£ are exponentially bounded. Note 

that since supp p T * p^ = supp [s/Q/kJF p^) = R and similarly for the 
second margin, we have supp p T = supp v T = R so that informational 
equivalence is guaranteed. In particular, if the moments of the mar- 
ginal distributions satisfy (j3J), then we know that the corresponding 
observables are informationally equivalent with position and momen- 
tum. Thus, the moment inequality (j3J) may be viewed as an operational 
sufficient condition for informational equivalence. 

Consider again the first marginal probability measure /i T * p^. We 
can now calculate the kth moment as 

(^ T *P?) [k]=j^( k \ T [k-n]$[n]. 

From this expression the moments of the position distribution can be 
solved recursively giving 

p?[*1 = if * p?) [*] - E - -]p?N- 

n=0 ^ ^ 

In other words, we are able to express the moments of the position 
distribution in terms of the operationally meaningful moments of the 
marginal statistics. If we know a priori that p^ is also exponentially 
bounded so that it is determined by its moments, then we have deter- 
mined uniquely the position distribution. This is the case, for instance, 
when p is a finite mixture of vector states from V, which is a dense set. 
A similar treatment can of course be carried out for the second margin. 

It is worth noting that if the measures p T and v T are exponentially 
bounded, then p T * p^ and z/ T * p£ are exponentially bounded whenever 
p^ and p^ are such. Therefore, the exponential boundedness of p T and 
v T guarantees that this method can be used for a dense set of states. 
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5. Conclusions 



We have shown that there is no direct connection between the in- 
formational completeness of a phase space observable and the state 
distinction properties of its Cartesian margins. More precisely, we 
have shown that it is not possible to infer informational complete- 
ness from the condition that the margins are informationally equiva- 
lent with sharp position and momentum observables. This shows that 
it is possible to determine the position and momentum distributions 
from the statistics of a single measurement even though the state is 
not uniquely determined. We have also demonstrated the converse fact 
that informational completeness does not guarantee the possibility of 
reconstructing position and momentum from the mere marginal sta- 
tistics, but one may occasionally need take the detour via full state 
reconstruction. 
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